A New Upper Bound on 2D Online Bin Packing 



Xin Han^, Francis Y.L. Chin^, Hing-Fung Ting^, Guochuan Zhang^, Yong Zhang^ 
^ Graduate School of Information and Technology, University of Tokyo 
hanxin . mail@gmail . com 
^ Department of Computer Science, The University of Hong Kong, Hong Kong 
{ chin, hfting, yzhang}@cs.hku.hk 
^ Department of Mathematics, Zhejiang University, China 
zgc@zju.edu.cn 

June 2, 2009 



Abstract 

The 2D Online Bin Packing is a fundamental problem in Computer Science and the determination 
of its asymptotic competitive ratio has attracted great research attention. In a long series of papers, 
the lower bound of this ratio has been improved from 1.808, 1.856 to 1.907 and its upper bound 
reduced from 3.25, 3.0625, 2.8596, 2.7834 to 2.66013. In this paper, we rewrite the upper bound 
record to 2.5545. Our idea for the improvement is as follows. 

In SODA 2002 [53], Seiden and van Stee proposed an elegant algorithm called H (g) B, comprised 
of the Harmonic algorithm H and the Improved Harmonic algorithm B, for the two-dimensional 
online bin packing problem and proved that the algorithm has an asymptotic competitive ratio of at 
most 2.66013. Since the best known online algorithm for one-dimensional bin packing is the Super 
Harmonic algorithm [5S], a natural question to ask is: could a better upper bound be achieved by 
using the Super Harmonic algorithm instead of the Improved Harmonic algorithm? However, as 
mentioned in |24| . the previous analysis framework does not work. In this paper, we give a positive 
answer for the above question. A new upper bound of 2.5545 is obtained for 2-dimensional online 
bin packing. The main idea is to develop new weighting functions for the Super Harmonic algorithm 
and propose new techniques to bound the total weight in a rectangular bin. 



1 Introduction 

In two-dimensional bin packing, each item {wi,hi) is a rectangle of width Wi < 1 and height hi < 1. 
Given a list of such rectangular items, one is asked to pack all of them into a minimum number of square 
bins of side length one so that their sides are parallel to the sides of the bin. Rotation is not allowed. 
The problem is clearly strongly NP-hard since it is a generalization of the one-dimensional bin packing 
problem [7] . In this paper we will consider the online version of two-dimensional bin packing, in which 
the items are released one by one and we must irrevocably pack the current item into a bin without any 
information on the next items. Before presenting the previous results and our work, we first review the 
standard measure for online bin packing algorithms. 



Asymptotic competitive ratio To evaluate an online algorithm for bin packing problems, we use 
the asymptotic competitive ratio defined as follows. Consider an online algorithm A. For any list L 
of items, let A{L) be the cost (number of bins used) incurred by algorithm A and let OPT[L) be the 
corresponding optimal value. Then the asymptotic competitive ratio for algorithm A is 

= lim supmax{A(L)/OPr(L)|OPT(L) = A;}. 



Previous work Bin packing has been well-studied. For the one-dimensional case, Johnson et al. |19j 
showed that the First Fit algorithm (FF) has an asymptotic competitive ratio of 1.7. Yao p8j improved 
algorithm FF with a better upper bound of 5/3. Lee and Lee [21] introduced the class of Harmonic 
algorithms, for which an asymptotic competitive ratio of 1.63597 was achieved. Ramanan et al. [23] 
further improved the upper bound to 1.61217. The best known upper bound so far is from the Super 
Harmonic algorithm by Seiden [25] whose asymptotic competitive ratio is at most 1.58889. As for the 
negative results, Yao [28] showed that no online algorithm has asymptotic competitive ratio less than 
1.5. Brown [1| and Liang [20] independently provided a better lower bound of 1.53635. The best known 
lower bound to date is 1.54014 [26] . 

As for two-dimensional online bin packing, a lower bound of 1.6 was given by Galambos [2]. The 
result was gradually improved to 1.808 [K], 1.857 [27] and 1.907 [1]. Coppersmith and Raghan [U] 
gave the first online algorithm with asymptotic competitive ratio 3.25. Csirik et al. |8j improved the 
upper bound to 3.0625. Csirik and van Vliet [lOj presented an algorithm for all d dimensions, where 
in particular for two dimensions, they obtained a ratio of at most 2.8596. Based on the techniques on 
the Improved Harmonic, Han et.al [T^ improved the upper bound to 2.7834. The best known online 
algorithm to date is the one called A ^ B presented by Seiden and van Stee [23], where A and B 
stand for two one-dimensional online bin packing algorithms. Basically A and B are applied to one 
dimension of the items with rounding sizes. In this seminal paper Seiden and van Stee proved that the 
asymptotic competitive ratio of -ff ® -B is at most 2.66013, where H is the Harmonic algorithm ^21] 
and B is an instance of the improved Harmonic algorithm. It has been open since then to improve the 
upper bound. A natural idea is to use an instance of the Super Harmonic algorithm [25] instead of the 
improved Harmonic algorithm. However, as mentioned in paper [24j, in that case, the previous analysis 
framework cannot be extended to Super Harmonic. 

We also briefly overview the offline results on two-dimensional bin packing. Chung ct al [6] showed 
an approximation algorithm with an asymptotic performance ratio of 2.125. Caprara [5j improved the 
upper bound to 1.69103. Very recently Bansal et al. [2j derived a randomized algorithm with asymptotic 
performance ratio of at most 1.525. As for the negative results, Bansal et al. [3] showed that the two- 
dimensional bin packing problem does not admit an asymptotic polynomial time approximation scheme. 

For the special case where items are squares, there is also a large number of results [9] [Ml [22l [TTl 
[121 [131 [18] . Especially for bounded space online algorithms, Epstein et al. [12] gave an optimal online 
algorithm. 

Our contributions There are two main contributions in this paper, 

• we revisit ID online bin packing algorithm: Super Harmonic, give new weighting functions for it, 
which are much simpler than the ones introduced in [25) . and the new weighting functions have 
interests in its own. 

• we generalize the previous analysis framework for 2D online bin packing algorithms used in [24j, 
and show that the new analysis framework are very useful in analyzing 2D or multi-dimensional 
online bin packing problems. 

By combining the new weighting functions with the new analysis framework, we design a new 2D online 
bin packing algorithm with a competitive ratio 2.5545, which improves the previous bound of 2.66013 
in SODA 2002 [24j. As mentioned in |24j, the old analysis framework does not work well with the 
old weighting functions in [25], i.e., the old approach does not guarantee an upper bound better than 
2.66013. This is testified in the following way: consider our algorithm, if we use old weighting functions 
with the old framework to analyze it, the competitive ratio is at least 3.04, and if we use the old weighting 
functions with the new framework, the competitive ratio is at least 2.79. 

Organization of Paper Section 2 will review the Harmonic and Super Harmonic algorithms as prelim- 
inaries. Section 3 defines the weighting functions for Super Harmonic. Section 4 describes and analyzes 
the two-dimensional online bin packing algorithm H SH+. Section 5 concludes. 



2 Preliminaries 



We first review two online algorithms for one-dimensional bin packing, Harmonic and Super Harmonic, 
wfiich are employed in designing online algorithms for two-dimensional bin packing. 

2.1 The Harmonic algorithm 

The Harmonic algorithm is a fundamental bin packing algorithm with a simple and nice structure, that 
was introduced by Lee and Lee [21] in 1985. The algorithm works as follows. Given a positive integer k, 
each item is immediately classified into one of k types according to its size upon its arrival. In particular, 
if an item has a size in interval (^^, \\ for some integer z, where 1 < i < k, then it is a type-z item; 
otherwise, it is of type-A:. The type-i item is then packed, using the simple Next Fit (NF) algorithm, 
into the open (not fully-packed) bin designated to contain type-i items exclusively; new bins are opened 
when necessary. At any time, there is at most one open bin for each type and any closed (fully-packed) 
bin for type-i is packed exactly with i items of type-z for 1 <i < k. 

For an item of size x, we define a weighting function Wh{x) for the Harmonic algorithm as follows: 



where H[L) is the number of bins used by the Harmonic algorithm for list L. 
2.2 The Super Harmonic algorithm 

The Super Harmonic algorithm |25j is a generalization of the Improved Harmonic algorithm and the 
Harmonic algorithm. Super Harmonic first classifies each item into one of A; -|- 1 types, where A: is a 
positive integer, and then assigns to the item a color of either blue or red. It allows items of up to two 
different types to share the same bin. In any one bin, all items of the same type have same color and 
items of different type have different colors. For items of type-i (i < k), the algorithm maintains two 
parameters Pi and ji to bound respectively the number of blue items and the number of red items in a 
bin. More details are given below. 

Classification into types Let ti = 1 > t2 > ■■■ > t/^ > t^+i = e > = be real numbers. An 
interval Jj is defined to be (tj+i, tj], for i = l,...,A;-|-l. An item with size x is of type-i if x G /j. 

Coloring red or blue Each type-i item is also assigned a color, either red or blue, for i < k. The 
algorithm uses two sets of counters, ei,...,ek and si,...,Sk, all of which are initially zero. The total 
number of type-i items is denoted by Sj, while the number of type-i red items is denoted by e^. For 
I < i < k, during the packing process, the fraction of type-i items that are red is maintained, i.e., 
ei = [ajSjJ, where ai, ...,ak € [0,1] are constants. 

Maximal number of blue items Let Pi = for 1 < i < A;, which is the maximal number of blue 
items of type-i which can be accepted in a single bin. 




The following lemma is directly from [21j. 



Lemma 1 For any list L, we have 



H{L) 



< 



J2WHip)+0{l) 



Space left for red items Let 6i = 1 — tiPi, which is the lower bound of the space left when a bin 
consists of /?j blue items of type-i. If possible, we want to use the space left for small red items. Note 
that in the algorithm, in order to simplify the analysis, instead of using Si, less space is used, namely 
D = {Aq, Ai, A;^}, as the spaces into which red items can be packed, where = Aq < Ai < • • • < 
Ak < 1/2 and K < k. Let A.^(i) be the space to be used to accommodate red items in a bin which 
holds Pi blue items of type-z, where function ^ is defined as {l,...,k} i— > {0,...,K} such that satisfies 
A^(j) < 5i. If = then no red items are accepted. 

For convenient use in our analysis in the next section, we introduce a function called ip{i), which 
gives the index of the smallest space in D into which a red item of type-? can be placed: 

ip{i) = mm{j\ti<Aj,l<j<K}. 

Maximal number of red items Now we define 7^. Let 7^ = if > Ak; otherwise 7^ = 
max{l, [Ai/ti\}, i.e., if Ai < ti < Ak, we set 7^ = 1, otherwise 7^ = [Ai/tjJ. 

Naming bins It is also convenient to name the bins by groups: 

{{i)\4>i = 0,l<i<k} 

{ii,'!Mi^O,l<t<k} 
{(?,j)|«i/0,l<j<fc} 
{{hj)\4>i / 0,aj 7^ 0,-fjtj < A<^(j),l < ij < k}. 

Group (i) consists of bins that hold only blue items of type-i. Group {i, j) consists of bins that contain 
blue items of type-i and red items of type-j. Blue group (i,?) and red group (?,i) are indeterminate 
bins currently containing only blue items of type-i or red items of type-j respectively. During packing, 
red items or blue items will be packed into indeterminate bins if necessary, i.e., indeterminate bins will 
be changed into 

The Super Harmonic algorithm is outlined below: 

Super Harmonic 

1. For each item p : z <— type of p, 

(a) \ii = k + 1 then use NF algorithm, 

(b) else Sj Sj + 1; if Cj < [ajSjJ then <— + 1; { color p red } 

i. If there is a bin in group (?, i) with fewer than 7^ type-i items, then place p in it. 

Else if, for any j, there is a bin in group (j, i) with fewer than 7^ type-i items then place 
p in it. 

ii. Else if there is some bin in group (j, ?) such that A^^^) > 7jfj, then pack p in it and 
change the bin into (j, i). 

iii. Otherwise, open a bin pack p in it. 

(c) else {color p blue}: 

i. if = then if there is a bin in group i with fewer than /3j items then pack p in it, else 
open a new group i bin, then pack p in it. 

ii. Else: 

A. if, for any j, there is a bin in group or (i,?) with fewer than /3j type-i items, 
then pack p in it. 

B. Else if there is a bin in group (?, j) such that A^(j) > ^jtj then pack p in it, and 
change the group of this bin into 

C. Otherwise, open a new bin (i, ?) and pack p in it. 



3 New Weighting Functions for Super Harmonic 



In this section, we develop new weighting functions for Super Harmonic that are simpler than the 
weighting system in p5]. The weighting functions will be useful in analyzing the proposed online 
algorithm as we shall see in the next section. 



3.1 Intuitions for defining weights 

Weighting functions are widely used in analyzing online bin packing problems. Roughly speaking, for 
an item, the value by one of weight functions is the fraction of a bin occupied by the item in the online 
algorithm. There is a constraint in defining weights for items for an online algorithm. Let K + 1 he the 
number of weighting functions. Let W^{p) be the weight of an item p, where 1 < i < K + 1. For any 
input L, the constraint is 

A{L) < ^^max^^ { ^ W\p)} + 0(1), (1) 

where A{L) is the number of bins used by algorithm A. 

Consider Super Harmonic algorithm. For 1 < i < A;, let li be the number of type-z pieces. For 
1 < i,s < A;, let -B(j ?), 5(7 be the number of bins in groups (i), (i, ?) and Then 

we have 

(1 - o:i)li 



E + E ^(M) + ^(.^) } = E + OH) 



and 



So, for each item with size x G /j, where i < A;, if we define its weight as below: 



1 - Qj _^ Oj 

Pi li ' 

then it is not difficult to see that the constraint ([1]) holds. However the above weighting function is not 
good enough, i.e., it always leads a competitive ratio at least 1.69103. 

The main reason is that for each bin in group (i, s) we account it twice, where 1 < i, s < k. Next 
we give some intuitions for improving the above weighting function. 

By ([2]) and ([3]), observe that 



and 



So, we have 



EEBis.)^E^ + oii). (5) 



v^v^D Ei Es + Ej Es -^(s.j) / ■s-^(l-ai)/i aj/j , 

2^2^%-) = 9 ^ 2^ 9/^.. +2^^ + 0(1 



Hence, for an item with size x € li, after packing, if there is a bin in group (i, s) and also a bin in group 
{s,i), then we can define its weight as below: 

1 - «i 

This is the main intuition to lead our weighting functions, which are given in the next subsection. 



3.2 New weighting functions 

Remember that in Super Harmonic, there is a set I? = {Aq, Ai, A^} representing the "free spaces" 
reserved for red items. Recall the two functions and ip{i) are related to free spaces and have the 
meanings as below: = j implies that free space Aj is reserved for red items in a bin consisting of /3j 
blue items of type-i, and ip{i) = j indicates that a red item of type-z could be packed in free space A>j. 

We are now ready to define new weighting functions. Items with size larger than e will be first 
considered. Let E be the number of indeterminate red group bins (?, i) when the whole packing is done. 

If £^ = 0, i.e., every red item is placed in a bin with one or more blue items, then we define the 
weighting function as: 

W\x)=^^^ ifxeli. (6) 

Pi 

Otherwise, E > implying that for some i, an indeterminate red group bin (?, i) exists after packing. 
Let e be the smallest red item in indeterminate red group bins. Assume r is the type of item e and 
j = (p{r). U2 < j < K then we define the corresponding weighting functions as follows: 



W^+^-^{x) 



( l-Qj _|_ OLi_ 

ft ^ "i-yi 



if X G /j (t){i) < J, and ip{i) < j 



ft 

l-a. 



li 



+ Si iixeh <P{i) > j, and (p{i) > j 
_|_ ^ if T (= 



(i) > j, and ip{i) < j 



If j = 1, we define 



W^+^{x) = < 







iixeli = 0, and ip{i) 
+ ^ iixeli (t){i) = 0, and ^{i) > 



ft 


7i 



if X G 7j (f){i) > 0, and ip{i) = 
a X e li (t){i) > and (p{i) > 



Note that in the above definitions, if 7^ = then we replace ^ with zero. For an item with size x G Ik+i, 
we always define W^{x) = for all j. 

Theorem 1 For any list L, we have 



where A{L) is the number of bins used by Super Harmonic for list L. 



Proof. Fix a list L. Let D be the sum of the sizes of the items of type-(A: + 1). By NEXT FIT, we know 
that the number of bins used for type-(/>; + 1) is at most D/{1 — e) + 1. 



Again, we use E to denote the number of indeterminate red group bins when all the packing is done. 
li E > 0, let e be the smallest red item in indeterminate red group bins. Assume r is the type of item 
e and j = ip{r). For 1 < z < A;, let /j be the number of type-i pieces. Let S(j^s), -B(j ?), 5(7 j) be the 
number of bins in groups (i), (i,s), (i, ?) and (?,i). 

To prove this theorem, we consider three cases. 

Case 1: If E = 0, i.e., Y2i ~ 0' every red item is packed in a bin with one or more blue items. 
Therefore we just need to count bins containing blue items: 

i s 

< y: w\x)+y:^^-^-p^+o{i) by® 

^ r -Pi 

W\x)+ W\x) + 0{1). 

Case 2: > 0, e is the smallest red item in indeterminate red group bins and its type is r and 
(p{r) = J > 2. Since every red item of type-z is placed in a final group bin (s, i), where ip{i) < j, we have 

<p{i)<j 

On the other hand, we have 

E = 0; (8) 

otherwise, e would have been placed into a bin (z, ?), where (f){i) > j. According to the Super Harmonic 
algorithm, for any type bin B^^^, we have 

c^ii) = 0. (9) 

Define 

^= E %^)' 

which is the total number of all the bins in groups (z, s) such that 0(z) > j and ip{s) < j. Then we have 

ML) < Y^ + E(% + %?)+%o) + EE%^) + o(i) 

i is 

= Y^ + E(^w+%7)+%.))+^+ E E%^)+ E Y^isA + oii) 

i <l){i)<j « •fi{i)>j * 

= E (% + %?)+E%^))+ E (%.) + E%o)+^ + o(i)- (10) 

The last inequality follows directly from ([T]), dS]) and ([9]). 
Then by the definition of variable X, we have 

E E%^)^^^dx< T.Bis,y 

j<<t>(i)<K s l<v{i)<j-l s 



Therefore, 

So, by (dni) and ([II]), we have 



< E (^w + %?) + E%^))+ E (%«) + E^(«.^) 



2 2 

(i)>j S (/3(j)<J s 



< 



D sr^ {^-Q.i)h sr^ Ojk sr^ (1 - ai)li ^-^ ajlj 

1 - e ^ A 7i 2/3j ^ 27j 

(/)(i)<j ■y3(i)>i 0(*)>i ¥'{«)<J 

D sr^ / (I - aj)/i Oi^jX ^p / (I - 0^)/^ Oj/j 



ip{i)<j v(i)>j 



+ E(^ + f)+E(^ + t) + oa) 

^(i)>j ^ ' *(i)>i ' 

¥'(*)>i 'P{i)<j 

W^+^~^{x)+ Yj W^-^^-^x) + 0{1) 

The second inequality follows directly from ([2]) and ([3]). 

Case 3. E' > and j = 1. The arguments are analogous with Case 2. According to the Super 
Harmonic algorithm, for any type of bin {i,s), we have ^p{s) > 1, where 1 < i, s < k and A; is a 
parameter defined in Super Harmonic. So, there is no such bin (z, s) with (p{s) < 1. Then we have 

ML) < E (% + %v) + E%^))+ E (%o + E^(^.o) + o(i) 

4>{i)<l s </3(i)>l s 

^ E E ^ + o{^) 

<^{j)=0 (/.(i)=0 " (/-(i)>0 

</3{j)=0 (/j(i)>0 ¥'(i)>0 

Therefore, we have A{L) < maxi<i<i^+i { ^^^^ W'{p)} + 0(1). □ 



4 Algorithm i7 <5'i7+ and Its Analysis 

In the section, we first review a class of online algorithms for two dimensional online bin packing, called 
H ® B [24j . Next we introduce a new instance of algorithm H ® SH+, where H is Harmonic and 
SH+ (Strange Harmonic+) is an instance of Super Harmonic. Then we propose some new techniques 
on how to bound the total weight in a single bin, which is crucial to obtaining a better asymptotic 



competitive ratio for the H ® B algorithm. Finally, we apply new weighting functions for SH+ to 
analyze the two-dimensional online bin packing algorithm H ® SH+ and show its competitive ratio 
at most 2.5545, which implies that the new weighting functions work very well with the generalized 
approach of bounding the total weight in a single bin. Note that as mentioned in [23] if we apply the 
weighting functions of SH+ derived from [25] directly to analyze algorithm H ^ SH+ then the upper 
bound cannot be improved. 



4.1 Algorithms H x B and H ^ B 

Now we review two-dimensional online bin packing algorithms H x B and H ® B [23], where H is 
Harmonic and B is Super Harmonic. 

Given an item p = (w, h), H x B operates as follows: 

1. Packing items into slices: If w > e then pack p into a slice of height 1 and width ti by H 
(Harmonic algorithm), where tj+i < w < ti] else pack it into a slice of height 1 and width e(l — 5y 
by H (Harmonic algorithm), where e(l — 5)*"*"^ < w < e(l — 5)* and (5 > is arbitrarily small. 

2. Packing slices into bins: When a new slice is required in the above step, we allocate it from a 
bin using algorithm B. 

H (S> B is a randomized algorithm, which operates as follows: before processing begins, we flip a 
fair coin. If the result is heads, then we run H x B; otherwise we run B x H, i.e., the roles of height 
and width are interchanged. Note that it is possible to de-randomize H ® B without increasing its 
performance ratio. For details, we refer to 



Theorem 2 If an online ID bin packing algorithm B has weighting functions W^{x) such that B{L) < 
maxi{J2xeL ^ij(^)} + 0(1). Then the cost by algorithm H ® B for input L is at most 

— ^ ( max { ^ T^|,,s(p)} + max { Vrij,^(p)}) + 0(1), 
and the asymptotic competitive ratio of algorithm H (x) B is at most 

— ^inax(max| J] t^j,(x)H^ij(y)| + max | WH{y)Wh{x)\] 

where 6 is a parameter defined in H B algorithm and X is a set of items which fit in a single bin. 



4.2 An instance of Super Harmonic SH+ 

As mentioned in [25j, it is a hard problem to find appropriate parameters in designing an instance of 
Super Harmonic, especially setting ti. The parameters in SH+ are found through a trial-and-error way 



and are defined as follows: 
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where ff{i) = 1.35(50 - i)/37{i - 12). 

Then we have seven weighting functions for SH+, i.e., 
1< f < 7. 



as defined in the last section. 



where 



4.3 Previous framework for calculating upper bounds 

In this subsection, we first introduce the previous framework for computing the upper bound of the 
competitive ratio of -ff ® SH+, then mention that the previous framework does not work well with the 
instance in the last subsection, i.e., the previous framework does not lead a better upper bound. 
Let p = {x, y) be an item. We define the following functions. 

Wh^nip) = WHix)Wh{y), Wh^nip) = WHiy)Whix), 

and 

^ WHix)W},iy) + Wiix)WHiy) 

Then we can obtain an upper bound on the competitive ratio R of algorithm H SH+ as follows 
by Theorems [T] and [21 where X is a set of items which fit in a single bin. 



R < TTT. max max 

- 2(1 -6) vx I i<i<7 



I E ^HxBiP)] + max I E ^BxHip)] 



< 



max 



(1 - 5) i<«,i<7,vx 



Y,W''\x,y)\ (12) 



The value of R can be estimated by the following approach. 

Definition 1 Let f be a function mapping from (0, 1] to 'P{f) is the mathematical program: 

maximize Exgx /(•^) subject to E^gx — ^' ^^^r all finite sets of real numbers X. We also use T'{f) to 
denote the value of this mathematical program. 

Lemma 2 \24^ Let f and g be functions mapping from (0, 1] to Let F = V{f) and G = V{g). 
Then the maximum of EpGX fi^{p))9{'^{p)) '^^^'^ finite multisets of items X which fit in a single 
bin is at most FG, where p is a rectangle and h{p) and w{p) are its height and width, respectively. 

In [21], / and g are defined as below: 
and 

g'-'[x)= sup — „. ., , ■ 
0<y<l f^'^y) 

By the above definitions, we have 

W''^ix,y)<f''^{y)g^'^ix), 

for ah < X < 1 and < y < 1. 

Remarks: As mentioned in [24], the old weighting functions [25] do not work well with the calcu- 
lating framework used in [251 124j . i.e., the previous framework with the old weighting function does not 
guarantee an upper bound better than 2.66013. 

4.4 A new framework for calculating upper bound 

In this subsection, we first generalize the previous analysis framework by introducing a new lemma and 
developing new functions for / and g in order to bound the total weight in a single bin. Then we apply 
our new weighting functions for Super Harmonic to algorithm H (8) SH+ and obtain a new upper bound 
for two-dimensional online bin packing. 

Lemma 3 maxyx | EpGX = maxyx | Epex where l<i,j <7 

Proof. By definition, observe that for any 1 < < 7, 

W'^^ix,y)=W^^%y,x). (13) 



Let X = {pi,P2, ■ ■ ■ ,Pm} be a set of rectangles which fit into a single bin, where pi = (xj, yj) is the i-th 
rectangle in X. If we exchange roles of x and y of pi to get new rectangles p'^ = {yi,Xi) for all i, then it 



is not difficult to see tfiat the new set X' = {pi^p^, ■ ■ ■ ,p'm.} ^^so a feasible pattern, i.e., all items can 
fit in a single bin. On the other hand, by equation we have 



where 1 < i,j < 7. There is a one-to-one mapping between X and X' in all the feasible patterns. 
Therefore, we have this lemma. □ 

New functions / and g: We define new functions / and g such that (i) Lemma [2] can be applied to 
bound the weight in a single bin, and (ii) the resultant bound is not too loose. The new functions / 
and g are defined as follows: 



f'^iy) = KjWHiy) + (1 - Kj)Whiy), 



where < Ajj < 1 and 



g^'-'ix) = sup 



0<y<l P'^{y) 



Note that in [2l|, Xij are 1/2 for all It is not difficult to see that the following inequality still holds 
although we have generalized the definition of the / function, 

W''^{x,y)<f'^{y)g''^{x) 

for ah < X < 1 and < y < 1. 

New approach of calculating V{f): In order to use Lemma [2] to obtain the upper bound on the 
competitive ratio R of algorithm H ® B, we need to calculate V{P'^) and 'P{g^'^). Let / be one of P'^ 
or (7*'-J for 1 < i, j < 7. In ^25j, Seiden wrote a programming to enumerate all the feasible patterns to get 
the bounds for V{f). Here, we give a simple approach by calling LP solver directly to estimate V{f), 
which can be modeled as the following mixed integer program (MIP): 



50 50 



max. 



/ = 'f^XiWi + {l-^Xi{ti+i + e)) ^ (1) 



=1 i=l 
50 

S.t. ^Xi{ti+i + e) <l, 

i=l 

Xj > 0, integer. 

where Xi is the number of type-i items in a feasible pattern, Wi is the weight for an item of type-i, which 
is decided by function /, i.e., Wi = f{p) if p € Since e > can be arbitrarily small, we cannot 

find an exact value for e. Therefore, we set e = and re-model the above MIP as follows. 

max. / = ^^l^x^w^ + {l-ZZlXit^+l)xT^ (2) 

s-t- Ylt^i XiU+i < 1, 

Xi<l, for 1 < i < 7, Xi<2, for 8 < i < 13, 

Xi < 3, for 14 < i < 15, Xi < i - 12, for 16 < i < 17, 

xis + xi9 < 6, Xi < i — 13, for 20 < i < 50, 

2x7 + ^15 ^ 3.9, 3x7 + 2xi3 + xi7 < 5.9, 

4x13 + 3x15 + 2:24 < 11.9, 5x7 + 3.53x11 + 1.47xi8 < 9, 

12x7 + 8x13 + 3x20 + 2:36 < 23, 9x7 + 6x13 + 2x2i + X30 < 17, 
Xj > 0, integer. 



Note that the new constraints do not ehminate any feasible solutions of MIP (1). For example, consider 
the constraint 5xi + 3.53xii + 1.47xi8 ^ 9. Since an item of type-7 has size larger than 0.5, an 
item of type-11 has size larger than 0.353 and an item of type-18 has size larger than 0.147, we have 
0.5x7 + 0.353x11 + 0.147x18 < 1. So, we have 5x7 + 3.53xii + 1.47xi8 < 10. It is not difficult to see 
that the following inequality 5x7 + 3.53xii + 1.47xi8 < 9 is equivalent to 5x7 + 3.53xii + 1.47xi8 < 10 
when X7, xii and xi8 are non-negative integers. For other constraints in MIP (2), the arguments are 
analogous. 

To solve MIP (2), we use a tool for solving linear and integer programs called GLPK |16j . We write 
a program to calculate W^'^ {x^y), g^'^{x) and f^'^{y) for each and then call API of GLPK to 

calculate V{P'^) and 'P{g^'^). The values of V{P'^) and V{g^'^) are shown in the tables in Appendix. 

Note that when we use Lemma[2]for the upper bound on the weight maxvx{l^pgx W^'^ {x^y)}, for 
all pairs (i, j), the calculations are independent. For different pairs (i, j), Ajj may be different. So, 
in order to get an upper bound near the true value of niaxvx{X]pex ^^'^ {^lU)}-: we have to select an 
appropriate Ajj. This can be done by a trial-and-error approach. 

Theorem 3 For all 5 > 0, the asymptotic competitive ratio of H ® B is at most 2.5545. 
Proof. According to the tables in Appendix, by Lemma and Lemma [21 we have 



max|EM^'''(p)| 
Sex J 



max 

vx 



W^^^ip) } < < 2.5539 



[■1,2 



maxj E = maxj E < < 2.5545. 

maxj E = max I ^ < r{f^'^)V{g^^^) < 2.5340. 

p£X ^ yX I J 

maxj E = max I ^ W^^\p)\ < V{f^''')V{g^^^) < 2.5364. 

pex ^ ^ p&x ^ 

For all the other by Lemma O we have 

maxj E < 'P{f''^)V{g''^) < r{f^'^)r{g^'^) < 2.5545. 



pex 



□ 



Remarks: If we use the weighting functions from [25j and the previous analysis framework, we 
find that the competitive ratio is at least 3.04. (run our programming 2DIISII.C like "./2DHSII-|-.exe 
old > yourfile") Even if we use the new weighting function, by the previous analysis framework, the 
competitive ratio is still at least 3.04, by running our programming 2DIISII.C like "./2DHSII-|-.exe newl 
> yourfile". We also find that if we use the old weighting function from [25] with the new analysis 
framework, the competitive ratio is at least 2.79. (run our programming 2DIISH.C like "./2DHSII-|-.exe 
old2 > yourfile") The reason is that: Lemma [2] does not work very well with the old weight function, 
i.e., the resulting value • G is away from the maximum weight of items in a single bin. 



5 Concluding Remarks 



When we use the tool for solving the mixed integer programs, there are two files which are necessary: 
one is the model file for the linear or integer program itself (refer to Appendix), and the other is the 
data file where the data is stored. We write a program to generate the data and then call the tool 
GLPK. (Actually we call API (Application Program Interface) of GLPK. To download the source file, 
go to: http:/ /sites. google. com/site/xinhan2009/IIome/files/2DHSH. c ). 

Our framework can be applied to 3D online bin packing to result in an algorithm H x H SH+ 
with its competitive ratio 2.5545 x 1.69103(f« 4.3198). 

Acknowledgments The authors wish to thank the referees for their useful comments on the earlier 
draft of the paper. Their suggestions have helped improve the presentation of the paper. 
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A Values of /*'^ and 





(1,1) 


(1,2) 


(1,3) 


(1,4) 


(1,5) 


(1,6) 


(1,7) 




0.500000 


0.500000 


0.540000 


0.550000 


0.565000 


0.565000 


0.600000 




1.598272 


1.598272 


1.605095 


1.606845 


1.609490 


1.609490 


1.615665 


rig''') 


1.598272 


1.597872 


1.574422 


1.581742 


1.585430 


1.587508 


1.575580 




2.554474 


2.553834 


2.527096 


2.541614 


2.551734 


2.555079 


2.545610 



ihj) = 


(2,1) 


(2,2) 


(2,3) 


(2,4) 


(2,5) 


(2,6) 


(2,7) 




0.500000 


0.500000 


0.530000 


0.550000 


0.565000 


0.565000 


0.600000 




1.597328 


1.597328 


1.597148 


1.597028 


1.596938 


1.596938 


1.596729 




1.609235 


1.598326 


1.586301 


1.595016 


1.602278 


1.604268 


1.589545 




2.570476 


2.553051 


2.533557 


2.547285 


2.558739 


2.561917 


2.538073 





(3,1) 


(3,2) 


(3,3) 


(3,4) 


(3,5) 


(3,6) 


(3,7) 




0.500000 


0.500000 


0.530000 


0.550000 


0.565000 


0.565000 


0.600000 




1.573676 


1.573676 


1.572837 


1.572777 


1.572732 


1.572732 


1.572627 




1.609235 


1.598326 


1.586301 


1.595016 


1.602278 


1.604268 


1.589545 




2.532414 


2.515247 


2.494992 


2.508604 


2.519954 


2.523084 


2.499762 





(4,1) 


(4,2) 


(4,3) 


(4,4) 


(4,5) 


(4,6) 


(4,7) 




0.500000 


0.500000 


0.535000 


0.550000 


0.565000 


0.565000 


0.600000 


np'') 


1.581245 


1.581245 


1.577140 


1.575380 


1.573621 


1.573621 


1.569515 




1.609235 


1.598326 


1.586855 


1.595016 


1.602278 


1.604268 


1.589545 


np'^w^n 


2.544594 


2.527344 


2.502692 


2.512755 


2.521378 


2.524510 


2.494814 



(iJ) = 


(5,1) 


(5,2) 


(5,3) 


(5,4) 


(5,5) 


(5,6) 


(5,7) 




0.500000 


0.500000 


0.535000 


0.550000 


0.565000 


0.565000 


0.600000 


V{P'^) 


1.585370 


1.585370 


1.580113 


1.577860 


1.575607 


1.575607 


1.570350 




1.609542 


1.598326 


1.587240 


1.595374 


1.602747 


1.604737 


1.589740 


V{P^^W'') 


2.551720 


2.533939 


2.508019 


2.517277 


2.525300 


2.528436 


2.496449 





(6,1) 


(6,2) 


(6,3) 


(6,4) 


(6,5) 


(6,6) 


(6,7) 




0.500000 


0.500000 


0.530000 


0.550000 


0.565000 


0.565000 


0.600000 


ViP'^) 


1.586853 


1.586853 


1.582237 


1.579160 


1.576853 


1.576853 


1.571468 




1.609785 


1.598326 


1.586682 


1.595657 


1.603117 


1.605107 


1.589894 


ViP'^Vig''^) 


2.554493 


2.536309 


2.510507 


2.519798 


2.527881 


2.531019 


2.498468 





(7,1) 


(7,2) 


(7,3) 


(7,4) 


(7,5) 


(7,6) 


(7,7) 




0.500000 


0.515000 


0.535000 


0.555000 


0.565000 


0.570000 


0.600000 


V{P'^) 


1.568686 


1.560602 


1.549821 


1.539044 


1.533655 


1.530958 


1.517143 




1.021572 


1.(309605 


1.602462 


1.612258 


1.G22822 


1.638219 


1.621359 


viP'^W'n 


2.543738 


2.511952 


2.483529 


2.481335 


2.488849 


2.508043 


2.464386 



B Model File for GLPK and Usage of Our Program 2DHSH+.C 

param I: =50; 



param c{i in l..I}>=0; 



param w{i in 1..I}; 

var x{i in 1..I}, integer, >=0; 

maximize f: sum{i in 1..I} w[i]*x[i] + (l-sum{i in 1..I} c[i]*x[i]) * 38/37; 

s.t. xO: sum{i in 1..I} c[i]*x[i] <= 1; 
xl : sum{i in 1..7} x[i] <= 1; 
x7: sum{i in 8.. 13} x[i] <= 2; 
xl4: x[14] <= 3; 
xl5: x[15] <= 3; 

xl6: x[16] <= 4; 
xl7: x[17] <= 5; 
xlB: x[18] + x[19] <= 6; 
y715: 2*x[7] + x[15] <= 3.9; 
y71317: 3*x[7] + 2*x[13] + x[17] <= 5.9; 
yl31524: 4*x[13] + 3*x[15] + x[24] <= 11.9; 
y71118: 5*x[7] + 3.53*x[ll] + 1.47 *x[18] <= 9; 
y7132036: 12*x [7] +8*x [13] + 3*x[20] + x[36] <=23; 
y7132130: 9*x[7] + 6*x[13] + 2*x[21] + x[30] <=17; 
others{i in 20.. 50}: x[i] <= i -13; 

end; 



Whene the parameters in Super Harmonic such as Oi, Pi, 7^ and (p{i) and ip{i) are given, we can 
calculate the weighting functions of Super Harmonic W^(-). Then the weighting functions 
for algorithm H ^ SH+ can be calculated as well as P'^{y) and g'''^{x). For each we call API of 

GLPK to solve V{f'^) and r{g''^). 

To use our program under linux system: 

• InstaU GLPK, 

• Compile: "gcc -o 2DHSH+.exe 2DHSH+.C -Iglpk" 

• Run: "./2DHSH+.exe new2 > yourfile" 

If there is an error message like "Could not load *.so" when you compile the source, then try to 
set "LD_LIRARY_PATH" as follows: "LD_LIRARY_PATH= $LD_LIRARY_PATH:/usr/local/lib", then 
"export LD_LIRARY_PATH". 



